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<D ; Abstract 

We derive four-dimensional relativistic three-body equations for the case of 
a field theory with a three-point interaction vertex. These equations describe 
■,—{ . the coupled 2 — > 2, 2 — > 3, and 3 — > 3 processes, and provide the means of 

calculating the kernel of the 2^2 Bethe-Salpeter equation. Our equations differ 
from all previous formulations in two essential ways. Firstly, we have overcome 
the overcounting problems inherent in earlier works. Secondly, we have retained 
all possible two-body forces when one particle is a spectator. In this respect, we 
show how it is necessary to also retain certain three-body forces as these can give 
rise to (previously overlooked) two-body forces when used in a 2 — > 3 process. The 
revealing of such hidden two-body forces gives rise to a further novel feature of 
our equations, namely, to the appearance of a number of subtraction terms. In 
the case of the ttNN system, for example, the NN potential involves a subtraction 
term where two pions, exchanged between the nucleons, interact with each other 
J> ! through the it-it t-matrix. The necessity of an input tt-tt interaction is surprising 

and contrasts markedly with the corresponding three-dimensional description of 
the ttNN system where no such interaction explicitly appears. This illustrates the 
somewhat unexpected result that the four-dimensional equations differ from the 
three-dimensional ones even at the operator level. 
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I. INTRODUCTION 

Relativistic quantum field theory is the underlying basis of essentially all models 
in both nuclear and particle physics. Yet our knowledge of how to solve any but 
the most simple of models using such a field theory is extremely limited. In this 
respect, it is useful to recall that, for the case of one- and two-particle systems, 
the way to sum all possible perturbation graphs is known, though only formally. 
In particular, the coupled problems of single particle dressing, three-point vertex 
dressing, and 2^2 scattering, is described in terms of the Dyson-Schwinger and 
Bethe-Salpeter equations. It is important to recognize, however, that despite the 

1 Accepted for publication in Nuclear Physics A 



1 



"text-book" nature of these equations, they do not constitute a solution of the 
one- and two-particle problem. The crucial point is that the essential ingredient 
in these equations is the kernel of the 2^2 Bethe-Salpeter equation, but this 
kernel is not known. Although one can calculate contributions to this kernel up to 
some order in the coupling constant, this is not a satisfactory method for strong 
interactions where it is necessary to have a non-perturbative approach. In this 
paper we develop integral equations for the coupled 2 — > 2, 2 — ► 3, and 3^3 
processes, which, in the case of negligible three-body forces, also provide exact 
closed equations for the unknown Bethe-Salpeter kernel. 

More generally, our goal is to formulate the field theoretic three-body problem 
where one-, two-, and three-particle states are explicitly coupled, and where three- 
body forces are taken as a given input (a reasonable first approximation is to 
take them to be zero). Unfortunately, the formulation of the four-dimensional 
three-body problem is not as straightforward as that for the two-body case. The 
basic problem is that " irreducibility" , the essential concept used to derive both the 
Dyson-Schwinger and Bethe-Salpeter equations, cannot be applied to the case of 
three-particles directly - as we shall see, quite a sophisticated additional analysis 
of diagram structure is needed. In fact, we find that in all the previous works on 
this subject, such analyses of diagram structure are either incorrect, incomplete, 
or missing altogether. Thus, despite the fundamental nature of the problem, there 
does not appear to be any satisfactory formulation of the three-body problem in 
four- dimensional relativistic field theory. The goal of this paper, therefore, is to 
provide just such a formulation. 

In order to develop equations for the field theoretic three-body system, we also 
utilise the concept of irreducibility. A Feynman diagram is said to be r-particle 
reducible if one can draw a continuous curve intersecting exactly r lines, at least 
one of which is an internal line, and each line being crossed just once, so as to 
separate intial from final states. A diagram that is not r-particle reducible, is said 
to be r-particle irreducible. The Bethe-Salpeter equation then follows from the 
following argument. Firstly, all 2 — ► 2 Feynman diagrams are classified according 
to their 2-particle irreducibility. The sum of all 2-particle irreducible diagrams is 
defined to be the potential for the process. All the rest of the diagrams have at 
least one 2-particle cut; moreover, it is clear from the topology of such diagrams, 
that there must be a unique left-most (or right-most) 2-particle cut. To the left 
of this left-most cut is a 2-particle irreducible potential, and the Bethe-Salpeter 
equation follows immediately. 

In a pioneering work, Taylor applied such a procedure to the three-body 
system. He recognised that, unlike the case for 2-particle cuts, 3-particle cuts may 
not be unique. Particular attention was paid to the case illustrated in Fig. |l](a). 
Here is shown a diagram that is 2- and 3-particle reducible, and where it is possible 
to cut three internal lines in two different ways, neither of which precedes the other. 
It was recognised that to eliminate such non-unique internal 3-particle cuts, one 
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Figure 1: Diagrams where there is no unique right-most 3-particle cut. (a) 2- 
particle reducible example, (b) 2-particle irreducible example. 



should remove the two-particle cuts first. Formalizing these considerations into 
his "last cut" lemma, Taylor used the scheme where 3-particle cuts are exposed 
only in 2-irreducible amplitudes, to derive four- dimensional relativistic equations 
for the field theoretic three-body problem. 

However, no special attention was given to the case illustrated in Fig. [l|(b). 
Shown is a diagram for the 2^3 process that is 3-particle reducible, but 2- 
particle irreducible. Shown also are two 3-particle cuts, neither of which is to 
the right of the other, i.e. there is no unique right-most 3-particle cut, despite 
the fact that two-particle cuts have been removed. This lack of uniqueness in the 
cutting procedure, necessarily involving at least one external particle, does not 
seem to have been taken into account by Taylor, or in a number of subsequent 
works that were based on his "last cut" lemma [2-5]; as a result, the three-body 
equations developed in all the above mentioned works suffer from the problem of 
overcounting. 

The only previous work where this type of non-uniqueness problem was ad- 
dressed, is the one by Tucciarone 0. As a result, he developed a description of the 
four- dimensional three-body problem where there is no overcounting of diagrams. 
However, his equations were not expressed in terms of subsystem two-body ampli- 
tudes, rather, an explicit perturbation expansion for the inhomogeneous terms is 
needed. This makes his equations unsatisfactory from a practical point of view. 

It should be noted, that the non-uniqueness difficulty does not arise for 3 — > 3 
processes where coupling to two-particle states does not take place, as for example 
in formulations of the relativistic four- dimensional three-nucleon problem 0. This 
problem also does not arise in time-ordered perturbation theory since, for example, 
the two cuts in Fig. |(b) would correspond to the two possible time orderings of 
the right-most vertices, both of which need to be taken into account. 

In the present work, we restrict ourselves to the case of a 3 field theory. This 
is the simplest case having wide applicability; however, the strategy we shall adopt 
is also applicable to other forms of the interaction. In particular, we overcome 
the non-uniqueness problem by a special procedure where, in cases like Fig. [l|(b), 
one of the two rightmost vertices is "pulled out" further to the right. Our equa- 
tions therefore do not suffer from the overcounting problem inherent in the above 
mentioned works. 



Figure 2: (a) A three-body force in the 3 — ^ 3 process, (b) The same three-body 
force, after particle 1 has absorbed particle 3, reveals a two-body process. 



A second, but equally important goal of this paper, is to take into account as 
many Feynman diagrams as possible, consistent with pair-like interactions. That 
this is not an altogether obvious task, is illustrated in Fig. 0. Fig. 0(a) shows a 
connected three-body force for the 3^3 process. Yet if we join any two initial 
or final legs into a single leg, i.e. if there is absorption through a 3 vertex, then 
the resulting diagram can be represented through a pair-like interaction, as shown 
in Fig. 0(b). Since the equations we seek couple 3 — > 3 to 3 «-> 2 processes, Fig. 
illustrates how careful one must be in neglecting three-body forces, otherwise one 
can inadvertently neglect important two-body contributions as well. This subtlety 
does not appear to have been noticed before, as all the References [1-5] have 
neglected three-body forces without comment. Note, that in an otherwise correct 
formulation, the neglect of such three-body forces, leading to the undercounting 
of pair-like interactions, would result in equations that have the wrong three- 
dimensional limit; thus, they would have neither the correct non-relativistic limit, 
nor would they satisfy three-body unitarity. 

In this work we present three-body equations that suffer neither from the over- 
counting nor the undercounting problems discussed above. A detailed derivation 
is presented for the case of three identical Bosons. The case of identical particles 
was chosen because it has all the complications of the most general case. It is also 
the most difficult case that can arise in practice, where, in particular, pair-like 
interactions are less straightforward to reveal. 

We also present the equations for the case of the ttNN system where only two 
of the particles are identical. This system has commanded much attention in the 
literature although mainly in the three-dimensional sector ||; in this respect, we 
note that the recently derived convolution equations of Ref. can be viewed as 
the three-dimensional limit of the four- dimensional equations we present in this 
paper. 

A novel feature of our equations is the presence of explicit subtraction terms 
necessary to avoid overcounting. In the ttNN equations, for example, one of the 
subtraction terms is the two-pion exchange NN potential where the exchanged 
pions interact with each other. As the corresponding three-dimensional equations 
have only the irN and NN interactions as input, it is remarkable that a consistent 
four- dimensional description of the irNN system also demands, as input, the ir-ir 
interaction. 



II. DERIVATION 

For the m — > n process the momentum space Green's function is defined by 
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where the S-matrix S is given in terms of the interaction Lagrangian Ci by 



S = Texp 



'J d A x Ci(x) 



(1) 



(2) 



For clarity of presentation, we shall usually replace momenta by their labels in 
an obvious way; thus we write Q(V . . . n', 1 . . . m) instead of G(p[ ■ ■ -p'mPi ■ ■ - Pm)- 
For momentum conservation 5 functions, we use a convention where, for example, 
8(1'2'3', 12) denotes the function 5(p[ + p' 2 + p' 3 — p\ — p-i). For the case n — m, 
it will also be useful to define Qo(l' . . . n' , 1 . . . n) as that part of Q{V ... n', 1 ... n) 
corresponding graphically to n completely disconnected lines connecting initial and 
final states. 

For specific one-, two- and three-particle Green's functions we shall write re- 
spectively 



<7(1', 1) = 0(1', 1) 
D(l / 2',12) = 0(1'2',12) 
G(l'2'3', 123) = ^(r'2'3',123). 



(3) 
(4) 
(5) 



We shall also need corresponding disconnected propagators where the momentum 
conserving 5 functions have been factored out; in these cases we write 



^(l , ,l)=d(l)5(l , ,l) 

A)(1'2', 12) = ^(1)^(2)5(1', 1)5(2', 2) 

Go(l , 2 , 3', 123) = ^(1)^(2)4(3)5(1', 1)5(2', 2)5(3', 3). 

For non-identical particles, we have that 



A)(l'2',12) =^ (r2',12) 
Go(l'2'3', 123) = £ (1'2'3', 123). 



(6) 
(7) 
(8) 



(9) 
(10) 



We shall work also with t-matrix amplitudes defined in terms of the above Green's 
functions by the relation . 

n rn 

T(l' . . . n', 1 . . . m) = J] d-,\j')g'{l' . . . n', 1 . . . m) J] d~\i) (11) 

j=l i=l 
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where Q' = Q if m ^ n, and Q' = Q — Go if m = n. Just as for propagators, we 
shall use specific symbols for the amplitudes of particular processes. Namely, we 
use symbol / for the dressed 1 — > 2 vertex, t for the 2-^2 process, T for the 1 — > 3 
process, M for the 2^3 process, and T for the 3^3 process. Thus 

/(1'2',1) = T(1'2',1) (12) 

t(l'2', 12) = T(1'2', 12) (13) 

r(l'2'3', 1) = T(l'2 / 3', 1) (14) 

M(l'2'3 / , 12) = T(l'2'3', 12) (15) 

T(l'2'3', 123) = T(l'2'3', 123). (16) 

The basic three-point vertex, without a momentum conserving S function is defined 
as h, thus 

/(1'2',1) = M1'2',1)5(1'2',1). (17) 
For the disconnected 2 — > 3 process it will be sufficient to define 

F^l'2'3',12) = /(l'3',l)d 2 - 1 (2)5(2 / ,2) (18) 
F 2 (l'2'3',12) = /(2'3',2)dr 1 (l)5(l',l). (19) 

To save on notation further, we shall usually suppress particle labels completely. 
Thus we denote any quantity A(V . . . n', 1 . . . m) simply by A. Note that the order 
of labels in this definition is important. As occasionally a different order of labels 
will be needed, it will be useful to introduce the left and right "particle exchange" 
operators Ljj and i?^; when acting on any quantity A, these operators exchange 
the ith and jth labels on the left and right set of external legs, respectively. For 
example, we have that 

Fi = Fi(l'2'3', 12), (20) 
R 12 F, = F^l'2'3',21), (21) 
L 13 F! = Fi(3'2T,12). (22) 

For reversed processes we use a bar notation; thus, / is the three-particle vertex 
for absorption, and M is the amplitude for the 3 — > 2 process. 

We shall also write integrals in symbolic form. For example, for any two quan- 
tities B and A, describing processes m — > k and k — > n respectively, we define AB 
by 

AB(p[ . . . p' n , Pl . . . Pm ) = J d'p'i . . . d A pl A(p[ ...p' n ,p'(... pl)B(p>{ . . . pi Pl ...p m ). 

(23) 

Thus an equation of the form 

C(p[...p' n ,p l ...p m )= J d A p'{ . . . d'p" k A{p' x ...p' n , p'{ . . . pl)B(p'{ . . . pl Pl ... p m ) 

(24) 
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will be written symbolically as just C = AB. 



A. Identical particle system 

We first consider the case of identical particles since this contains all the com- 
plications of the most general case. To be definite, we assume that the identical 
particles are Bosons and that the underlying interaction Lagrangian is given by 

C T (x) = ^(x). (25) 

Here A is the coupling constant and the factor 1/3! is chosen to cancel factors 
arising from the permutation symmetry of the fields in contractions when using 
Wick's theorem. Having defined the model in which we shall work, it will now be 
assumed that all three-point vertices and particle propagators are dressed, so we 
consider the skeleton Feynman diagrams of a given process. 

For identical Bosons, the Green's function Q of Eq. (|l|) is symmetric under the 
interchange of either initial or final labels. In this case we still retain all the above 
definitions, so the full propagators D, G, and t-matrices /, T, t, M, and T are 
likewise symmetric in their initial and final labels. Note, however, that the free 
propagators D and Go do not have a specific symmetry, and the function is 
symmetric under the interchange of % and 3 momentum labels, but does not have 
a specific symmetry under the exchange of labels 1 and 2. 

Quantities without a specific symmetry, will sometimes need to be symmetrized 
explicitly. Thus, for example, we have that 

^ J D (l'2 / ,12) = ^o(l , 2 , ,12) (26) 
p 

]T G (l , 2 , 3', 123) = 0o(l'2'3', 123) (27) 
p 

where the sum is over all permutations of either the initial or final momenta. 
More generally, we shall indicate sums over permutations of initial momenta by 
the letter R (right), and sums over permutations of final momenta by the letter L 
(left); however, we use the letter P whenever it makes no difference which sum, R 
or L, is taken. Quantities symmetrized in one of these ways are indicated by the 
appropriate superscript. Thus, for example, if A has three initial and final legs, 
then 

A R (1'2'3', 123) = A(l'2'3', 123) + A(l'2 , 3', 213) + A(l'2'3', 321) + . . . (28) 
A L (1'2'3', 123) = A(l'2'3', 123) + A(2T3', 123) + A(3'2'l', 123) + . . . (29) 
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with similar expressions holding for A having any number of legs. In general, we 
can write 



A R 

A L 
A p 



R 

E^ 

L 

J2A = A R = A L . 



In the case of three-particle states, the symbols R c , L c , and P c will be used to 
indicate sums over the corresponding cyclic permutations of particles. 

We shall utilise the idea of irreducibility described in the Introduction. The 
maximal irreducibility of an amplitude will be indicated by a superscript in round 
brackets; thus by we mean the sum of all possible 3 — > 3 skeleton Feynman 
graphs that are simultaneously 1-, 2-, . . . , r-particle irreducible. 

The symmetry property of identical particle amplitudes necessitates the inclu- 
sion of various counting factors, as well as sums over permutations of particle labels 
whenever amplitudes do not already have the required symmetry. In the Appendix, 
we discuss the origin of such counting factors, and give particular examples of how 
to symmetrize amplitudes. As a result, we shall utilise such counting factors and 
symmetry operations in the discussion below, but with only a minimal comment. 



1. The 3 — > 3 process 

We begin by examining the structure of the full amplitude T for the 3^3 
process. We can expose the one-particle cut through the equation 



T = T (1) + r (1) #r (1) 



(30) 



which is illustrated diagrammatically in Fig. |[ In the same way, we can expose 
the two-particle cut in 



T (D = T (2) +M (2) [ l jDo+ 1^(1)1^^(2) 

= T^ + ^M^[D^ + D t^D ]M^. 



(31) 



1 





f — s 

1 









Figure 3: Graphical representation of the equation T = + T^-'gT^ 1 '. 



Note that the factors of 1/2 in the first of the above two expressions are associated 
with the symmetric nature of the amplitudes and under the interchange 
of the exposed two particles; more generally, exposing n particles between appro- 
priately symmetric amplitudes will result in a counting factor of 1/n! - see the 
Appendix for details. By exposing the first vertex in amplitude we can write 
it as 

T« = ±M®D f (32) 

which is illustrated in Fig. |j. Note that because all our vertices are dressed from 
the beginning, the 2 — > 3 amplitude in Eq. (|3~2"D is 1 and 2-particle irreducible. 
Substituting Eq. (|31|) into Eq. fl3~0|) and using Eq. (p2) we obtain 



(33) 



T = T m + -M^DM^. 
4 



where we have used the fact that 



t = t 



(i) 



fgf 



and 



D 



D£ + D tD . 



(34) 



(35) 



In Eq. (pop we can treat t in two ways. In order to compare with some earlier 
works we would treat i as a given input to the equations. In the second way we 
would expose three-particle states in t and thereby obtain coupled equations for 
the 3 — > 3 and 2 — > 2 processes. Let us assume for the moment that t is known. 
Then the only unknowns in Eq. (|33|) are the amplitudes and . 

For the amplitude T^> any three-particle cut is unique as there are no two- 
particles states in this amplitude. This means that we can immediately write the 
three-particle Bethe-Salpeter equation for this amplitude 



T (2) =K + ±-kg q tw 

3! 



(36) 



1 




Figure 4: Graphical representation of the equation = ^M^D f 
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where K is the sum of all 3 — > 3 diagrams that are simultaneously 1-, 2- and 
3-particle irreducible. We write the disconnected part of K, indicated by subscript 
d, in terms of 1- and 2-particle irreducible two-body potentials v (= t^) 

K d {l'2'3', 123) = J2 V C2'3', 23)^ (1)6(1', 1) (37) 

L C R C 

where L c and R c indicate that the sum is taken over cyclic permutations of the left 
labels (1'2'3') and right labels (123) respectively (note that the sum is restricted to 
cyclic permutations because the potentials v are already symmetric in their labels). 
Note that 

t {l) =v + \vD t {l) . (38) 



Defining by 



we have that 



V 1 (l'2'3 , ,123) = w(l'3 / ,13)^ 1 (2)(5(2 / ,2) (39) 
V 2 (l'2'3', 123) = v{2'3',23)d- 1 {l)5{l',l) (40) 
V 3 (l'2'3', 123) = v{l'2',12)d- 1 {3)5{3',3), (41) 



Kd = £(K + ^ + V 3 ) (42) 

Pc 



where it makes no difference over which labels, left or right, the cyclic permutations 
are taken. The subscript i of Vi can be considered as labelling the interacting pair. 
In this respect our convention throughout this paper shall be as in the above 
equations, i.e., 1 labels the pair (13), 2 labels the pair (23), and 3 labels the pair 
(12). Note that the definitions of Fi given earlier in Eqs. ( |l8l) are consistent with 
this convention. Further, we have that 

K = J2(Vi + V 2 + V 3 )+K c (43) 

Pc 

where K c is the connected part of K, and as such, is a three-body force. Defining 

V = l -{V 1 + V 2 + V 3 )+ l -K c (44) 
2 o 

so that 

tf = £v; (45) 
p 

Eq. (|36|) can be written as 

T {2) = Y^V + VG T {2) . (46) 
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In obtaining Eq. ( pE6|) we made use of the fact that is symmetric; it is therefore 
gratifying to see that this derived equation now implies the symmetry of T^ 2 ' 
explicitly. 

In this work, we would like to be able to neglect three-body forces, but as 
discussed in the Introduction, this needs to be done very carefully In Eq. ([46]), we 
can in fact neglect K c in the expression for V; however, as we shall see shortly, this 
can be done only because there are no two-body states in the amplitude T^ 2 \ As it 
stands, Eq. flip]) does not have a compact kernel so a Faddeev rearrangement needs 
to be performed. To facilitate this, one can recast Eq. (|46|) into a more familiar 
form (for the scattering of three non-identical particles) by defining the t-matrix 
T through the equation 

f = V + VG f (47) 

where 

T®=Y,T. (48) 
p 

Eq. (|47D can then be transformed to give equations of standard Faddeev form. Al- 
though the introduction of T is the familiar way of handling Eq. (|46|), we show that 
in fact one can work with Eq. fl46|) directly. To do this, it is more straightforward 
to work with Green's functions. We thus write 

= qp + GqT^Gq, (49) 
= G% + G VG {2 \ (50) 

Both these equations for G^ and Eq. (|46"D for differ from the corresponding 
ones for non-identical particles only in the symmetrization of the inhomogeneous 
term; moreover, in the equations for G^ 2 \ this symmetrization is of the free Green's 
function Go and is thus particularly simple. Without the K c term, Eq. (|50|) is given 
by 

G m = G p + l Go ( Vl + y 2 + vr 8 ) G ra (si) 

We define the components Tj by 

TiG = ±ViG® (52) 



so that 



and 



G® = G^ + Y l G T i G . (53) 



T (2) = ^ + T2 + ^ ( 54) 

Using Eq. (|53| ) in Eq. (|52D , we obtain 



TiGq — —ViGq + —Vi ^2 GoTjGo, (55) 
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(1 - Iv-GoW = ^V^G, 1 + -ViGo J2 Tj, (56) 
and thus the Faddeev equations for the components T iy 

Ti = \$ + \uG Q Y. T s (57) 

where 

U = Vi+ ^ViGoU (58) 

Note that the inhomogeneous term of Eq. ( |57D involves a sum over right (but not 
left!) permutations. Comparing with Eq. fl38D , we see that the tj are explicitly 
given by 

ti(l'2'3',123) = t (1) (l'3 , ,13)rf- 1 (2)5(2',2) (59) 
t 2 (l'2'3',123) = t (1) (2 , 3',23)rf- 1 (l)5(l',l) (60) 
t 3 (l'2'3',123) = t (1) (l'2 , ,12)rf- 1 (3)(5(3',3). (61) 



2. The 2 — > 3 process 

Although Eqs. fl57|) and (54) determine the 1- and 2-particle irreducible am- 
plitude T^ 2 \ we still need to determine the amplitude before the full 3 — > 3 
amplitude T can be specified by Eq. (|33|) . 

Now to obtain the full connected 2 — » 3 amplitude M c , we write it in terms of 
one-particle irreducible and reducible parts: 

M c = M { p + Y [1) gf. (62) 

Similarly we can write 

MP = MW + -MWD Q t®. (63) 
Using Eqs. Q32|) and (|34]), we then obtain 

M c = M c (2) + M (2) ^D t 

= MP~D t + MP(l + ~D t), (64) 

where the disconnected part of M*- 2 ** is clearly given by 

Mf=£(F 1 + F 2 ). (65) 
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It is seen, therefore, that all we need now is the connected amplitude Mj?\ and 
then both the coupled 3 — > 3 and 2 — > 3 processes will be fully determined. 

In the case of amplitude M^ 2 \ the very first possible three-particle cut is not 
unique since the amplitude contains two-particles in the initial state. However, 
here this non-uniqueness is simpler to handle than in the case of amplitudes with 
two-body intermediate states, since, apart from the two initial particles, Mj?' 
has no further two-particle states. In particular, we can avoid the overcounting 
problems associated with the non- uniqueness of the first three-particle cut, by 
simply "pulling out" a vertex / on one of the two initial particles. Furthermore, the 
clearest way to see the structure of is to temporarily expose the bare vertex 
function f° that is buried inside this chosen vertex / of the skeleton diagram; 
namely, we write 

MP = \[TMG F?\c (66) 

where F® is defined as in Eq. ( fT%| ) but with f° replacing /. Eq. ( Efh is illustrated 
in Fig. [|. At this stage it is more straightforward to work in terms of Green's 
functions, so we use Eq. (f49|) to write the alternative form of Eq. (|6SD, 



= \[G^F«] C . (67) 

We also introduce the Green's functions Gi defined by the equations 

d = D^{l'3',13)g{2',2) (68) 

G 2 = DW(2'3', 23)0(1', 1) (69) 

G 3 = £>W(1'2',1%(3',3). (70) 

Note that the Gi are symmetric under the interchange of labels corresponding to 
the interacting pair. It is easy to see that 

G t = GP + ^GoViGi (71) 

= G? + G tiG (72) 





7 — ^ 








Mo 





Figure 5: Graphical representation of = \[T^G Q F^\ 
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where the symbol Pi indicates a sum over permutations of (left or right) momenta 
of the corresponding interacting pair: 

Gq 1 = (1 + P 13 )G = P P (1'3', 13)2(2', 2) (73) 
Gq 2 ee (1 + P 23 )G = P P (2'3', 23)^(1', 1) (74) 
Gq 3 ee (1 + P 12 )G = P P (1'2', 1%(3', 3). (75) 

It is then evident that 

Co Pi = \g.FI (76) 
and using Eq. ( ffll) to write (1 — \GqVi)Gi = G Pl , we obtain 

(1 - ^G V 1 )G Q F 1 = = G F? (77) 

from which follows that 

Pi° = (1 - ^iG )Pi. (78) 



Using this result in Eq. (|6 

^^K^fl-^iW (79) 
Writing the alternative form of Eq. (J50[) , 

G( 2 ) = G p + G (2) UG , (80) 
the previous expression becomes 

= \Gv 1 GM{V 2 /2 + V 3 /2 + K C /6)G F 1 . (81) 

where the term Gq 1 GqFi has been dropped as it has no connected part. To obtain 
a scattering equation for M^> we multiply Eq. fl8"T|) on the left by 1 — VGq and 
note that Eq. (EDI) implies 



(G 1 - 10G< 2 ) = G^Cft = £ 5(1', 1)5(2', 2)5(3', 3). (82) 

p 

Thus 

MP = VG MP + ~ £04/2 + V 3 /2 + P c /6)G Pi. (83) 

1 L 

With the derivation of this equation, we have, at least formally, solved the double 
counting problem. Eq. ( [33|) is an integral equation whose solution provides 
and consequently, the full Green's functions for both the 2 — > 3 and 3 — > 3 process; 
of course, this equation has kernels V{ which are non-compact, but this is an easily 
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Figure 6: Two views of the same Feynman graph, (a) Viewed this way, the graph 
belongs to K c GqFi where K c is a three-body force, (b) Viewed this way, the graph 
belongs to V1G0F2 where V\ is a two-body potential. 



overcome formality. In a three-dimensional time-ordered approach, we would get, 
instead of Eq. (j83|) , the equation 

MP = VG M^ + W( V 2/^+y3/^+K c /Q)G F 1 + W( v i/ 2 + v ^/ 2 + K c/Q)G F^ 

1 L L 

(84) 

which, in contrast to Eq. (83), is symmetric with respect to the particle labels 
1 and 2. The seeming asymmetry of Eq. (|83|) implies that there must be two- 
body forces, hidden within the three-body force K c , that restore the symmetry. 
Revealing these hidden two-body forces is therefore of great conceptual importance, 
quite apart from questions of their importance in practical calculations. 

The amplitude K c is one-, two- and three-particle irreducible, and it is con- 
nected. In this sense, it is a genuine three-body force and it would be tempting to 
follow the common procedure of neglecting this term. This, however, would be a 
major error: as illustrated in Fig. |6], the combination K c GqFx, entering Eq. ( j33|) , 
in fact contains pair-like interactions when viewed as having the very first vertex 
as F 2 . 

Let us specify more precisely how this dual nature of the one diagram in Fig. 
^ manifests itself in our theory. To do this, we need to relate our four- dimensional 
three-body problem to the corresponding formulation in three-dimensional time- 
ordered perturbation theory. This can be done as follows. Let J- (pi . . . p n , q± . . . q m ) 
be some Feynman diagram contributing to the Green's function Q of Eq. ([]]), and 
define J-(E, pi . . . p n , qi . . . q m ) by the expression 

„ n 

^(£, Pl ...p n , qi ...q m ) = J 5(E-Y,P-)dp° 1 ...d P n 

8=1 
m 

T{j>x ...p„,qi... q m )5(E - £ g°) dq* . . . dq° m . (85) 

i=l 

Note that the integrations in Eq. flS5[) are in fact with respect to the relative 
energies of the Feynman diagram. It can then be shown, that T is the sum of 
all possible orderings of a time-ordered perturbation graph topologically identical 
to J 7 , together with additional terms due to the advanced part of the two-time 
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Green's function [ ID| . Thus if we take the set of Feynman diagrams specified by 
K C G F U then 

KcG Q Fi = Js{e-pI-pI- p° 3 )dp Mdp° 3 

K c G F l (p lP2 p 3 , qi q 2 )5(E - g° - q°)dq[dq° (86) 

includes the sum of all possible orderings of corresponding time-ordered pertur- 
bation theory graphs; amongst these are diagrams whose graphical representation 
is given by Figs. ||(a) and (b). Now in time-ordered perturbation theory, the di- 
agram of Fig. |(b), is unambiguously a two-body rescattering in the presence of 
a spectator, and is therefore included in any three-dimensional approach; indeed, 
neglecting this diagram while keeping other two-body contributions would lead to 
violation of three-body unitarity. For these reasons it is of utmost importance that 
what we have called pair-like contributions in K c GqFi be retained. 

To examine more closely which pair-like diagrams of the form ViG^F 2 (2 = 1,3) 
also belong to the three-body force term K c GqFi, we expand Vi by order of the 
interaction: 

V = V} 2] + V} A] + . . . (87) 

where V} n ^ is the ra-th order potential term. The example of Fig. |6] shows already 
that the fourth-order term, v} A \ gives a diagram that is a member of K c GqF\. 
Clearly, therefore, all higher order terms of Vi will likewise give diagrams belonging 

[21 

to K c GqFi. This leaves only the lowest order term, VI , as the single exception as 

[21 [21 

it does not belong to K c GqF x . The corresponding terms V{ G F 2 and V 3 G F 2 , 
illustrated in Fig. [Tj, are just appropriately symmetrized versions of the basic graph 

5(1'2 , 3', 12) = 5(1'2'3', 12)d(l")d(2") h(l'l", l)h(-l" - 2", -2 , )/i(3'2", 2), (88) 

shown as the bottom left diagram of Fig. [?]. 



2' , 2 2 

V^GqF 2 = 3/ 



1' 



+ 1' 
1 3' 



V 3 2 GqF 2 





to 




1" 



2 3' 

+ 1' 
1 2' 



Figure 7: Terms of the form V.GoF 2 that do not belong to K c GqFi. The bottom 
left diagram is denoted by B in the text. 
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Taking the above discussion into account, we can write 

K C G F 1 = [(Vi - V} 2] ) + (V 3 - V 3 [2] )} G F 2 + K lc G F x (89) 

where K\ c is made up of graphs which do not give rise to a two-body interaction 
when multiplied from the right by GqF\ with a subsequent "pulling out" of ver- 
tex F 2 . Note that the combination Kx c GqFi, despite its notation, is completely 
symmetric; moreover, diagrams belonging to Ki c GqFi do not give rise to pair-like 
interactions when either F\ or F 2 is pulled out to the right. This suggests a possi- 
ble classification of 2 — > 3 diagrams: (i) those that give pair-like interactions when 
either of the two vertices is pulled out - the only possible case is shown in Fig. |7|, 
(ii) those that give pair-like interactions when just one of the two vertices is pulled 
out; this case is given by 



{Vx - vn + (V 3 - V 3 [2] )\ G F 2 + [(V 2 - VD + (V 3 - V^)\ G F 1 , 
and (iii) those that do not give rise to pair-like interactions - Ki c GqF\. Defining 

Bo = ^(v} 2] + V 3 [2] )G F 2 , (90) 

and using Eq. ([89]), we can write Eq. ( |S3"D as 

(1-VG )MP = l^[\(V2 + V 3 )G F 1 + ^(V 1 + V 3 )G F 2 -B ] + ^K lc G F 1 (91) 

where, because of the left permutation sum, B can also be replaced by IB. By 
symbolic form, this result differs from the corresponding result in time-ordered 
perturbation theory, just in the subtraction of the term B . This can be understood 
because, in our covariant case, the right most vertices are not time-ordered relative 
to each other, so that Bq is contained in both the first and second terms on the rhs 
of Eq. (0). Thus one subtraction of Bq is necessary to avoid overcounting. The 
potential V on the lhs of these equations still contains the three-body force K c as 
given by Eq. (|44]) . Now, however, we can safely neglect this term as it is multiplied 
on the right by a connected amplitude, and thus cannot give rise to pair-like 
interactions. With the understanding that all genuine three-body forces are to be 
neglected, Eq. ( |9~TD provides a linear integral equation for the amplitude where 
all potentials are pair-like. We now need to perform a Faddeev rearrangement in 
order to obtain a compact kernel. This we cannot do directly to Eq. (^T|) because 
of the sum over left permutations. 



r[2U . ,,[2] 



By analogy with Eq. (47), one way to proceed would be to define the amplitude 
M satisfying the equation 

(1 - VG )M = i[i(F 2 + V 3 )G F 1 + hy x + V 3 )G F 2 - 2B] (92) 
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Figure 8: Graphical representation of the sums (a) El ^GqFx = 2V 2 G Fi + 
2R l2 V l G F 2 + 2V 3 G F 1 and (b) El V 3 G F 1 = 2V 3 G F 1 + 2V 2 G F 1 + 2R 12 V 1 G F 2 ; 
only the cyclic permutations of the left hand labels are shown. 



where 

MP=Y,M, (93) 

L 

so that Eq. fl9"2| ) is the same as Eq. ( ]9~T| ) but without the sum over permutations. 
We may then express M in terms of the set of equations 

M = Y.^-B, (94) 

i 

Mi = ViGo^M + Pt], (95) 

where 

Pi = \ft.\ fo = \Fi\ fo = \(F 1 + F 2 ). (96) 
A Faddeev rearrangement then gives the equations 

M t = - ±B) + hiGo^ Mj- (97) 

This method corresponds to a reduction to the different particle case; however, as 
we now proceed to show, this reduction is not necessary as the sum over permuta- 
tions in Eq. ( |9~ID can also be explicitly evaluated. 

The evaluation of El V2G0F1 and El V 3 G Fi in Eq. (|9"T|)is most easily seen 
diagrammatically, as illustrated in Fig. It is evident that 

EWi = EWi- (98) 

L L 

In a similar way, we can consider the graphical representation of the terms El V1G0F2 
and El V^GqF^. In this way, we deduce that 

£[(V 2 + V 3 )G F 1 + (V, + V 3 )G F 2 ] 

L 

= AV 1 G F 2 R + 4V 2 G F* + 4V 3 G F* (99) 
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where, on the rhs, the potentials Vi appear as in the different particle case, without 
left-hand side permutation sums. 

Substituting the last equation into Eq. fl9"T]) , we obtain 

= VyGaF* + V 2 G F« + V 3 G F« - B L + ^K lc G Q F x . (100) 

For the description of the 2 — > 3 process, we can drop the terms involving K c 
and Ki c in this equation. It must be emphasized, however, that once these three- 
body force terms are neglected, the resulting amplitude should not be used 
to calculate the 2 —>■ 2 amplitude directly, i.e. through a multiplication on the left 
by FiGq] such a procedure would certainly miss important pair-like interactions 
coming from the term FiGqK c and, as we shall see shortly, also from the term 
FiGqK 1c GqFi. At this stage we consider only the 2 — > 3 process so the neglect 
of K c and K lc is justified. To lowest order, Eq. (|100|) gives the expected result 
= B L , thus providing a simple check of our result so far, and illustrating the 
necessity of retaining the subtraction term B L . 

We can now perform the standard Faddeev rearrangement of Eq. ( |100| ) by 
defining 

M c (2) = Y. M *- bL i ( 101 ) 

i 

M t = KGo^AfW + Oi] (102) 

where 

ai = F<f ; a 2 = Ff; a 3 = F? or F 2 R . (103) 

With the help of Eq. (|5BD, Eq. ( |102| ) can be written in terms of two-body t-matrices 
as 

M % = UGoiai -B l ) + \Y j UGqMj. (104) 

—i ■ i ■ 



3. The 2 — > 2 process 

In the above discussion, the 2 — > 2 amplitude t was assumed to be a known 
input. In the following discussion, t will not be assumed to be known, rather, we 
use the given field theoretic model to write an integral equation for t. 

The first steps, separating out the one- and two-particle cuts in t, have already 
been performed in Eqs. ( [34]) and (|38|). We therefore begin by examining the struc- 
ture of the 2 — » 2 potential v. Since v is connected and one-particle irreducible, 
both of the final state particles in v have separate left-most vertices. In order to 
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expose three-particle states, we shall "pull out" the left-most vertex on the first 
particle. This is just the same procedure as was used in Eq. (|66|) for exposing 
three-body states in M^ 2 \ We recall, that in the case of we pulled out the 

bare vertex fo- However, we soon showed that this resulted in an expression, Eq. 
(|79|), where it is the dressed vertex / that has effectively been pulled out. This 
result has its origin in the general property of quantum field theory where a careful 
summation of Feynman diagrams leads to an effective perturbation theory where 
only skeleton diagrams need be considered. 

To underscore this equivalence of working with bare and dressed vertices, we 
shall this time pull out the left-most dressed vertex on the first particle, thus we 
write 

v = ^FiG M', (105) 

which is illustrated in Fig. Here M' consists of all the diagrams of M (2) except 
those that lead to dressing of vertices and particle propagators in the above ex- 
pression for v (these restrictions are also sufficient to guarantee the connectedness 
of v). Examples of such diagrams belonging to M&> but excluded from M' are 
given in Fig. [H| 

Writing M' in terms of its connected and disconnected parts, 

v = ±F 1 G M' c + ±F 1 G M' d (106) 

one can simply use Eqs. (|100|) and flS5|) for and M d respectively, keeping 
only those terms that do not give rise to propagator or vertex dressing; in this 



sense, we write Eq. ( |106|) as 

v= l -{F lGo Mi 2> } ND + l{PiGM 2) } ND (107) 

(2) 

where subscript ND stands for "no dressing". The part of v resulting from M d 
is thus given by 

v° PE = \ {hG,Mf} ND = F lGo F« (108) 





7 — ^ 


M' 


/ 







Figure 9: Graphical representation of v = \FiG$M' . 
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(a) (b) (c) 

Figure 10: Example of diagrams, belonging to M^ 2 \ that upon multiplication 
from the left by F X G , lead to dressing graphs. Diagram (a) leads to propagator 
dressing, (b) and (c) lead to vertex dressings. 



which is just the properly symmetrized one-particle exchange (OPE) potential. 
The full potential v is thus 



v = v OPE 



+ ^-hG^GoF* + V3G0F* - X> + \Ki c G q F x ) 

+ \f 1 G q {\v 2 + 1 -V, + \k c )G M^) . (109) 
2 2 2 b ) nd 

To display the pair- interactions resulting from K c , we use the analogue of Eq. 

F X G Q K C = £ F 2 G, [(Vt - V} 2] ) + (V 3 - V 3 [2] )] + F 1 G K' lc . (110) 

Rc 

Note that K[ c ^ K lc . Moreover, F x GqK' 1c G MP does not have pair-like rescat- 
tering contributions (with only one particle as spectator) as can be easily seen by 
using the lowest order contribution for M^> . Similarly, we define the analogues of 
Eqs. ( |88D and fl90|) for the 3^2 process by 

B{1'2\ 123) = 5{l'2\ 123)d(l")d(2") h{l'l", l)h{-l" - 2", 2)/i(2'2", 3). (Ill) 

and 

B = ±F 2 G (V} 2] + V 3 m ) (112) 
respectively. Using Eqs. ( 110 ) and ( |112j ) in Eq. ( |109| ), we obtain 

v = v OPE + l - |f 1 G (^ 2 G'oF/ ? + V 3 G F* " E 5 + \Ki c G g F x ) 

+^F X G Q {V 2 + V 3 ) + \f 2 G {V x + V 3 ) -B + \f x G K' 1c \G M?A . (113) 

2 2 O J ND 

At this stage, Mj 2 ) in this expression is still defined by the exact expression of Eq. 
( |100| ). However, in contrast to the case of Eq. ( |107| ), the factors multiplying 
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in Eq. ( |113| ) ensure that the three-body force contributions to do not give rise 
to hidden two-body contributions in Eq. ( |113| ). We may therefore safely neglect 
such three-body contributions and thus have Mj 2 -* described by Eq. ( |101|) rather 
than Eq. (|100|). 

In Eq. ( M terms of the form ViG M^ may be eliminated using Eqs. ( |102| ) 
and (|103| ), and the symmetry of may be used to reduce B to 2B. Thus, 
neglecting "safe" three-body forces, we obtain 

v = v OPE + X - j-AGo + \f 1 G K 1c G F 1 + F 1 G M 2 + F 2 G M 1 

+(F 1 + F 2 )G M 3 - 2BG MP - F 2 G V 1 G F 2 R - F 2 G V 3 G F*} ^ . (114) 

We still need to identify pair-like interactions that might be present in the term 
involving the three-body force Ki c ; that is, we seek hidden expressions of the form 
FiGoVkGoFj, including cases where the initial and final momenta are permuted. 
In this respect, we note that all such pair-like interactions can be expressed as one 
of the four terms FiG V 3 G Fj = 1,2); for example, 

FiG V 2 G Fi = FiG V 3 G Fi, 
F 2 G ViG F 2 = F 2 G V 3 G F 2 , 
Li 2 FiGqV 3 GqFi = F 2 GqV 3 GoFi, 
etc. 

(of course cases like F 2 G V 2 G F 1 do not arise as our vertices are dressed from 
the beginning). We note, therefore, that the last two terms of Eq. ( |114| ) may be 
simplified: 

~(F 2 G V 1 G Ff t + F 2 G V 3 G F{ t ) = F 2 G V 3 G (F 2 + F x ) (115) 

= F 2 G V 3 G F 2 R (116) 



the result being illustrated in Fig. [Ll| Note that this term is not symmetric under 
the interchange of its left labels. This implies that there must be a compensat- 
ing term in Eq. ( |1 14|) that restores the symmetry. Indeed we will see that this 
compensating term can be found in FiGoKi c GqF\. 

Central to the term FiGqKi c GqFi is the three-body force K\ c . Although K\ c 
has been defined [by Eq. QSPp] so that no graph of Ki c G Fi is expressible as vertex 
F 2 followed by a two-body potential, the full term F^qKicGqFi, nevertheless, 



involves just such a two-body potential. In Fig. |12]we give an example of a graph, 
belonging to Ki c , which, on being sandwiched between FiGq and GoFi, is of the 
form F 2 GoV 3 GoF 2 . It is important that such pair-like interactions are not neglected 
in the theory. 
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Figure 11: Graphical representation of F2GoVaGoF 2 R . 



To identify all such pair-like interactions in FiGqKi c GqFi, it is sufficient to 
examine the four terms FiGoV^GoFj = 1,2). In this respect, we note that 
terms of the form V 2 G F 1 or V 3 G Fi cannot belong to K 1c GqFi since K lc is a 
three-body force. Likewise, terms of the form V1G0F2 or V3G0F2 cannot belong to 
K\ c GqFi just by the property of K\ c given below Eq. (]89f) . Consequently, out of the 
above mentioned four terms, two of them, FiGqV$GqFi and F1G0V3G0F2, cannot 
be of the type FiGqKi c G$Fi. This leaves only two candidates that are consistent 
with the definition of K\ c , namely F2G0V3G0F2 and F2G0 V3G0F1. However, even 
these have contributions which do not belong to FiGqK]_ c GqFi. In particular, we 
need to exclude from F 2 G V 3 G F2 and F 2 G V 3 G Fi terms of the form FiGqV^GqFi 
and FiG V 3 G F 2 , as well as terms corresponding to vertex dressing (that these are 
the only terms to be excluded may be checked explicitly by using Eq. ( |114j ) to 
specify V3). The only such terms in our theory are given by the amplitudes W, X 
and Y defined by Fig. [13], together with the same amplitudes but with the labels 
on the two external right legs interchanged. The diagram of Fig. |TB|(a), defining 
amplitude W, corresponds to two-particle exchange where the exchanged particles 



(b) 



Figure 12: (a) Example of a graph belonging to the three-body force K\ c . (b) The 
graph of (a) sandwiched between operators F\Gq and GqFi, the resulting graph 
belongs to FiG Ki c G F 1 , but is nevertheless of the two-body type F2G0V3G0F2. 
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(a) (b) (c) 

Figure 13: Graphs that need to be subtracted from the two-body rescattering term 
F2G0V3G0F2 as they do not belong to FiGqKi c GqFi. (a) Interacting two-particle 
exchange, W; in the s-channel, the circle represents the amplitude t' = — v OPE . 
(b) Crossed two-particle exchange , X. (c) Vertex dressing graph, Y. 



interact through the interaction t', where 

t r = t (l)_ v OPE (n7) 

Here the subtraction of v OPE is needed to eliminate contributions to vertex dress- 
ing. The diagram of Fig. 0(b), defining amplitude X, corresponds to a "crossed 
two-particle exchange" potential; here the exchanged particles do not interact. The 
diagram of Fig. |l3|(c), defining amplitude Y, corresponds to vertex dressing and 
must be excluded from our final expressions. In summary, we can write that 

^F 1 G K lc G F 1 = F 2 G V 3 G F R -W R -2X-Y R + ^F 1 GoK 2c G Fi. (118) 

where K 2c is a three-body force which when sandwiched between F\Gq and GqF\, 
does not result in a pair-like interaction. Note that the amplitude X is symmetric 
and is thus equal to its exchange term; this and other symmetry aspects of this 
expression are discussed in detail in the Appendix. 

On substituting the results of Eqs. ( |116| ) and (|118| ) into Eq. ( |114| ), the unsym- 
metrical terms of Fig. [11] are eliminated, and we obtain 

1^ 



v ope + |_[ J p lGo M 2 + f 2 g oMi + (Fx + F 2 )G M 3 ] 



--F^qY.B- BGqMP - W R -2X -Y R \ (119) 
2 l J ND 

where we have neglected the term involving the three-body force K 2c . In this 
expression, we can further evaluate 

^FxGq Y,B = X+(1 + L l2 R 12 )Y. (120) 
1 L 
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As written, Eq. ( |119| ) specifies that dressing terms inside the curly brackets 
need to be removed. This we can now do explicitly by identifying all possible 
dressing terms, and then including subtraction terms in order to remove their 
contribution. Some vertex dressing contributions have already been revealed in 
the Y amplitudes of Eqs.( |119|) and ( |120| ). The only other vertex dressings arise 
from the first iteration of the terms F^G^Mj in Eq. ( |119| ). Thus, with the Mj 
evaluated to first order using Eqs. (|104|) , 

-[FiG M 2 + F 2 G M 1 + (Fi + F 2 )G M 3 ] ^ owest orc i er ) 

= -^[FiG^GoF^ + F 2 GotiG F 2 R + (Fx + F 2 )G t 3 G Fl i ] 
= (Fx + F 2 )G t 3 G (Fx + F 2 ) = J2 F 2 G t 3 G F 2 . (121) 

LR 

The vertex dressing terms arise from the one-particle exchange potential part of 
£3, that is, 

Y,F 2 G v° PE G F 2 = Y,Y, (122) 

LR LR 

and these must be subtracted from Eq. ( |119| ). Putting this together, the potential 
for the 2 — > 2 process is 

v = v OPE + \[F X G Q M 2 + F 2 G Mx + (Fx + F 2 )G M 3 ] - BG Q M® 

/Li 

-W R -3X -J2 Y - ( 123 ) 

LR, 



B. The 77 NN system 

Since the equations we have derived for identical particles contain all the com- 
plications of the most general case, it is now essentially a formality to write down 
the corresponding equations for other systems involving three-point vertices. Here 
we would like to apply our formalism to the case of the ttNN system which, be- 
cause of its fundamental nature, has been under intensive investigation for many 
years ||. It is thus assumed that we are dealing with a field theory of nucleons 
and pions where the interaction is given by a tcNN vertex. 

The discussion of the previous section does not need to be modified in any 
essential way. However, one does need to take into account that now one has only 
two identical particles, and moreover, that these two are Fermions. In this respect, 
we need to modify our definitions of the various symmetry operations introduced 
in the previous section. Although the symbols L and R still have the meaning of 
sums over permutations, now the sums extend only over the two nucleon labels, 
and with the additional condition that the "exchanged" term enters with a minus 
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sign. Choosing the convention that label 3 always refers to the pion, Eqs. and 
(|29f) are now modified to 

A R (1'2'3', 123) = A(l'2'3', 123) - i(l'2'3' 213) (124) 
A L (1'2'3',123) = A(l'2'3', 123) - A(2T3', 123) (125) 

with similar relations holding for quantities having only nucleons in initial or final 

states. 



1. ttNN -> nNN 

The amplitude T for the ttNN — > 7riV7V process is given, as for the identical 
particle case, by Eq. (|33|). The 2-particle irreducible amplitude T^ 2 ^ is again the 
sum of Faddeev components T& as in Eq. ([54]) ; however, the equations for these 
components are now given by 

Ti = t? + UG Y,Tk, (126) 
T 3 ^t 3 + h 3 G ^T k , (127) 

with 

U = % + ViGoU, (128) 
h = V 3 + ^V 3 G t 3 , (129) 

where, in the above, i — 1,2, and = 1,2,3. Note that £j is the 7rA^ amplitude 
where the pion interacts with the zth nucleon, the other nucleon being a spectator. 
Similarly t 3 is the NN amplitude in the presence of a spectator pion. 



2. NN -> vrAW 

With three-body forces neglected, the integral equation for the NN — ► 7rAW 
2-particle irreducible amplitude M^ 2 ) is given by a slightly modified form of Eq. 
fllUUp , namely 

- (Vx + V 2 + l -V 3 )G,M^ = VxGoF* + V 2 G F 1 R + V 3 G Q F? - B LR (130) 

where no 1/2 factors appear in front of V\ and V 2 since they correspond to potentials 
between distinguishable particles (irN). The amplitude B now takes on the form 
shown in Fig. [14|, and sums over both left and right nucleon exchanges are now 
also needed. 
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Figure 14: The term B in the case of the ttNN system. 



To effect the Faddeev rearrangement of Eq. ( |130| ) we define the amplitude 
components by appropriately modified Eqs. (|104|) 

Mi = ViG [MP + Fj*\, (j^i) (131) 
M 3 = V 3 Go[^MP + F 1 r ], (132) 
MP = Y, M k-B LR - (133) 

k 

In the above equations, as in the two equations below, we use the convention that 
i, j = 1,2 while k = 1, 2, 3. The rearranged equations are then given by 

M 3 = hG Q {F^- l -B LR )+ l -Y, t ^ M ^ (134) 
Mi = t i G (Fj t — B LR ) + ^ UGoM k (j^i). (135) 



3. NN -> NN 

For NN scattering, there can be no contributions that are 1-particle reducible 
in the specified model where all interactions are generated by the irNN vertex. 
The full NN t-matrix, tNN, is thus given by the Bethe-Salpeter equation of Eq. 
(|38|), which in the present case takes the form 

1 

tNN = vnn + 2 Vnn ^°^ nn - (136) 

Here -Do is the dressed two-nucleon propagator, and the NN potential, vnn, is 
specified in a similar way to Eq. ( p.23| ), namely 

vnn = v% E + F 1 G M 2 + F 2 G M 1 + (Fj + F 2 )G M 3 

-B LR G MP - W R - - W R N -X R - Y LR . (137) 

where v N ^ E is now the NN one-pion exchange potential and the terms W n7r , W^n, 
Wnni X r and Y are illustrated in Fig. 113. Each of the three different W amplitudes 
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(a) (b) 




• • •• • 

v ^ ^ i v : 1 



(d) (e) 
Figure 15: Graphs that need to be subtracted from the NN potential in the 
relativistic formulation of the irNN system, (a) W^, (b) W n N, (c) Wnn, (d) X, 
and (e) Y. The small circles represent the t-matrix amplitudes: (a) t^, (b) t!^ N , 
and (c) ijvjy — v^ E , each being the appropriate version of Eq. ( |117| ). 



involves the exchange of two particles that interact with each other through the 
appropriate version of the t-matrix t' of Eq. ( |117|) . Thus the pions of amplitude 
W n7T interact through the full ti-tt t-matrix t nn , the amplitude W w n involves 7r- 
A" scattering through the so-called background term, t^ N , and Wnn involves NN 
scattering through t/viv — v nn E '• Note that Wnn does contain u-channel one-pion 
exchange, although t-channel one-pion exchange is forbidden. 

It is important to note that the amplitude tNN enters our equations both as 
an input and as an output. In this respect, we note that the t-matrix £ 3 , defined 
formally by Eq. (|6lD , in our case is just the NN t-matrix tNN i n the presence of a 
spectator pion: 

i 3 (l'2'3', 123) = t NN (l'2', 12)c/; 1 (3)(5(3', 3). (138) 
Similarly, Eq. (|41|) becomes in our case 

y 3 (l'2'3', 123) = v N n(1'2', 12)C 1 (3)5(3', 3), (139) 

and Eq. (|129| ) can thus be identified with Eq. ( |136| ). Since £3 enters as an input 
to Eqs. (|127|) and (|134j) , so does t N N- Similarly, t N N is needed as an input to 
calculate Wnn- At the same time, t N N is the output amplitude of the Bethe- 
Salpeter equation Eq. ( |136| ). Ideally, the input and output tNN need to be self- 



consistent. This might be achievable through an iterative process where the output 
of one iteration becomes the input of the next iteration. Alternatively, a great 
numerical simplification is afforded by sacrificing the self-consistency and using an 
externally constructed tNN for the input. This may be reasonable if the observables 
of interest are not very sensitive to the details of the AW channel. 
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With the possible exception of £aw, the two-body scattering t-matrices entering 
the calculation of the W^-amplitudes form part of the externally constructed input 
to our equations. It seems remarkable that the 7T-7T t-matrix is a necessary input 
to our equations. No such input is needed for the corresponding three-dimensional 
theory of the irNN system, nor has it appeared in previous formulations of four- 
dimensional irNN equations. In the present work, the appearance of this input 
arises through the careful retention of all nN interactions where one nucleon is a 
spectator. 

SUMMARY 

We have derived four- dimensional relativistic three-body equations based on 
3 field theory. These equations couple the amplitudes for 2 — > 2, 2 — > 3, and 
3^3 processes, and can be considered as the natural extension of the Dyson- 
Schwinger and Bethe-Salpeter equations to the three-particle sector. With the 
assumption that three-body forces can be either neglected or given explicitly, these 
three-body equations form a closed set of equations for the above amplitudes, and 
provide the exact non-perturbative solution of the field theoretical problem. In 
particular, these equations provide an answer to the long standing question of 
how to determining the kernel of the 2 — > 2 Bethe-Salpeter equation in a non- 
perturbative way. If three-body forces are not given, then in the spirit of our 
approach, one could expose four-particle states in the three-body forces, and thus 
express them in terms of four-particle equations where now four-body forces are 
either neglected or explicitly given. Although one can, in this way, envisage a 
whole hierarchy of equations, involving ever higher multi-particle forces, there is a 
large amount of evidence that three-body forces are small for most processes; thus, 
the three-body equations developed in this paper, should form a solid theoretical 
starting point for the investigation of relativistic few-body processes. 

Two features distinguish our approach from previous attempts to formulate the 
relativistic four-dimensional equations. Firstly, we have overcome the overcounting 
problem, present in all previous works that use the so-called "last-cut lemma"; to 
do this, we used a procedure of explicitly pulling out one initial vertex in the 
2^3 amplitude. Secondly, we have retained all possible two-body interactions 
taking place in the presence of a third spectator particle. In previous works, it was 
apparently not realised that three-body forces in the 3 — » 3 amplitude can give rise 
to such two-body forces in the 2^3 amplitude. By retaining appropriate three- 
body forces, we avoid the serious problem of undercounting pair-like interactions 
which would arise if all three-body forces in the 3^3 amplitude were to be 
neglected. 

Another interesting aspect of our four- dimensional equations is that they dif- 
fer at the "operator level" from the corresponding three-dimensional three-body 
equations of time-ordered perturbation theory. This is in contrast to the case 
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for the two-particle system where the four- dimensional Bethe-Salpeter and three- 
dimensional Lippmann-Schwinger equations differ essentially in the dimension of 
integrals involved, but otherwise have the same operator form. For the three-body 
problem, this difference is due to the presence of a number of subtractions terms 
in our four- dimensional equations which are not present in the three-dimensional 
equations. These extra terms are necessary to avoid overcounting and arise as a 
direct consequence of including all the above-mentioned possible two-body interac- 
tions. We note that all the previous works that suffer from overcounting [1-5] have 
four- dimensional equations that are of the same operator form as the corresponding 
three-dimensional equations. 

It should be emphasised that our four- dimensional three-body equations are not 
only of theoretical interest. With the ever increasing power of current computers, 
these equations should also provide a timely practical tool for calculating few-body 
hadronic processes in nuclear and high-energy physics. 
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APPENDIX 

Here we would like to explain the origin of the counting factors appearing in 
expressions of Section II dealing with identical particles. 

The most common such factor is the inverse factorial 1/n! associated with n- 
particle intermediate states connecting two amplitudes which are symmetric with 
respect to label interchanges of these n particles. Although this factor can be 
derived by applying Wick's theorem to the Green's function of Eq. ([!]), its appear- 
ance is not a feature of a four-dimensional approach, rather, it is due to the basic 
symmetry property of identical particle amplitudes. Indeed, just the same factor 
appears in three-dimensional approaches where, in fact, the origin of this factor 
can be more easily seen. In a three-dimensional approach, the Green's function 
corresponding to the four- dimensional one of Eq. (P is given by 

G = < 0\a(p' n ) . . . <Pi) E+ \ H a\T>i) ■ ■ ■ a j (p m )\0 > (Al) 

where H is the Hamiltonian and the creation and annihilation operators obey the 
commutation relations (for Bosons) 

[a(p),aV)] =*(p-p'); [a\p) } a^p')] = [a(p), a(p')] = 0. (A2) 
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Just the same creation and annihilation operators appear in the expansion of the 
fields <p(x) in Eq. (p]) (here there is no subscript on since all particles are identical): 

<l>{x) = I , % Hp)e~^ + at(p)e**]. (A3) 

For three-dimensional Green's functions of Eq. (|A1|) , the states 

|pi . . . p n > = a f (pi) . . . a f (p n ) |0 > (A4) 

form a basis of the full Fock space of identical particle states. As is well known 
|TT| , the closure relation in this Fock space is given by 

00 f 1 

1° >< °l + zl / rf Pl • • • rf P™ — |Pl • • • Pn >< Pl • • • Pn\ = 1, (A5) 
n=l J n - 

and here one can explicitly see the origin of the 1/n! factor for n-particle interme- 
diate states. 

We note that together with the condition < 1 > = 1, Eqs. (|Aj) or Eq. 
iH) define the normalization of the states |pi...p n >. This normalization is 



particularly convenient for the case of disconnected graphs as we now demon- 
strate. Let Gi(oii, Pi) and Qiia'^P'^) denote two four-dimensional Green's func- 
tions that make up the two disconnected parts of the disconnected Green's func- 
tion Qd(&i(Xj, PiP'j)- Here we have used a short-hand notation for the momen- 
tum labels: e.g. a, = {ax, ■ • •} denote the left-hand momenta of Q\ and 
ctia'j = {«!, a 2 , • • • , a[, a' 2 , ■ ■ ■} denote the left-hand momenta of Note that 
Qi, Q2 and are symmetric functions of their labels. Now it is straightforward to 
show that 

g^a^PiP^)^ GifaPifa^ffj) (A6) 

where the sum is over distinct terms resulting from all possible exchanges of mo- 
mentum labels, as indicated. The above mentioned normalization ensures that this 
expression appears with no extra counting factors. The result of Eq. (|A6|) can be 
obviously generalized to disconnected Green's functions made up of any number 
of disconnected pieces, there again being just a sum over exchange of labels with 
no extra counting factors. In Section II, Eqs. (0), (pTf), ( p?|) and ( j65| ) can be 



considered as special cases of this result. It must be emphasized, that although 
Eq. (|A(J) is presented for four- dimensional Green's functions, the structure of the 
result, in terms of a simple sum over label exchanges, is due to the normalization 
and symmetry properties of the underlying Fock states of Eq. (A3). Thus the same 



equation holds for three-dimensional Green's functions, although in that case, the 
simple product of Q\ and Qi in Eq. (|A6|) needs to be replaced by a convolution 
integral fT2|j . 
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Figure 16: Illustration of how to determine the properly symmetrized amplitudes 
for (a) non-interacting two-particle exchange, and (b) interacting two-particle ex- 
change. 



We now would like to discuss the various counting factors entering the expres- 
sion of Eq. (|118|) . This will at the same time demonstrate how counting factors 
can be determined for most other expressions. 

The factor of 1/4 entering the left side of Eq. ( |118| ) is associated with the vertex 
functions F± and F\ as they expose two identical particles each. 

Turning to the r.h.s. of Eq. ( |118| ), let us show that the amplitude X of Fig. 
0(b), constructed simply from four symmetric vertices and without any extra 
counting factors, is just the properly symmetrized non-interacting two-particle 
exchange amplitude in our model. In Fig. [l6|(a) we demonstrate how such a sym- 
metrized two-particle exchange amplitude is constructed. The factor of 1/16 in 
Fig. |i~6l(a) results from factors of 1/2 for each of the four vertices bracketing the 
symmetrized four-particle propagator. The factor of 1/16 could also have been 
obtained by first symmetrizing the products FiF 2 and FiF 2 using Eq. ( |A6|) ; in 
each case there would be 6 terms, each giving an equal contribution. One would 
then need to multiply twice by 1/4!, since a symmetrized product of two vertices 
exposes 4 particles in intermediate state. Thus (6/4!) 2 = 1/16, as before. Now 
consider the sum over permutations P of the 4-particle propagator in Fig. |T6|(a): 
there will be 24 terms in this sum, however, not all of them will lead to a two- 
particle exchange amplitude. In fact it is easily seen that 8 of these permutations 
will result in disconnected dressing graphs, leaving only 24 — 8 = 16 that are proper 
two-particle exchanges. Now all of these 16 contributions are equal to each other, 
so taking into account the previous factor of 1/16, we obtain that the properly 
symmetrized non-interacting two-particle exchange amplitude is just given by X, 
the amplitude defined by the single diagram of Fig. 0(b). In a similar way, let's 
find the symmetrization factor necessary in front of the interacting two-particle 
exchange amplitude W defined in Fig. |13|(a). The properly symmetrized interact- 
ing two-particle exchange is constructed as in Fig. |l6|(b). We again have the factor 
of 1/16 due to the four vertices. This time, however, two of the four particles in 
the center of the diagram interact. In Fig. [I6|(b), the symmetrization sum over V 
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and R! is, as in Eq. flM|), over exchanges of left and right labels leading to distinct 
contributions. Now out of the 4! x 4! possible left and right permutations, the 
switching of labels 2 <-» 3, 2' <-> 3' and 1' <-> 4', leads to a repetition of already 
included terms, so they must be excluded. The sum over V and R' thus extends 
over only 4!4!/8 terms. However, from all these terms, there will be cases where 
the (interacting) legs 2 and 3 will both be connected to the same vertex, thus 
leading to a vertex dressing graph. Such cases must be excluded. Thus, effectively, 
there are only 4 possibilities for legs 2 and 3: one leg connected to either of the 2 
legs of the top left vertex, the other being connected to either of the 2 legs of the 
bottom left vertex. Similarly there are 4 possibilities for the 2' and 3' legs of the 
middle diagram to connect to the right-hand vertices. In total there can therefore 
be only 4 x 4 = 16 interacting two-particle exchanges, and taking into account the 
initial 1/16 factor, we see that the properly symmetrized interacting two-particle 
exchange is given by the amplitude W of Fig. |T3](a) with no counting factors. 

This kind of reasoning can be used for other types of diagrams involving iden- 
tical particles. 
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